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We investigate the eflects on cosmological density perturbations of dark radiation in a Randall- 
Sundrum 2 type brane-world. Dark radiation in the background is limited by observational con- 
straints to be a small fraction of the radiation energy density, but it has an interesting qualitative 
effect in the radiation era. On large scales, it serves to slightly suppress the radiation density per- 
turbations at late times, while boosting the perturbations in dark radiation. In a kinetic (stiff) era, 
the suppression is much stronger, and drives the density perturbations to zero. 



I. INTRODUCTION 
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' In brane-world cosmology, which has emerged in the context of recent developments in M theory, the observable 
universe is a l-|-3-dimensional "brane" surface embedded in a l-l-S-l-d-dimensional "bulk" spacetime. Fields and 
particles in the non-gravitational sector are confined to the brane, while gravity propagates in the bulk. Simple 
phenomenological brane-world models are the cosmological generalizations [1-4] of the Randall- Sundrum 2 model [5], 
which has a self-gravitating brane with Minkowski geometry embedded in an infinite part of 5-dimensional Anti de 
Sitter space (AdSs). Gravity on the brane is prevented from "leaking" into the infinite extra dimension at low energies 
—{ , via "warping" of the metric by the negative bulk cosmological constant, A5 — — 6/^^, where i is the curvature radius 
K*" ' of AdSs. Modes of the 5D graviton have an effective mass on the brane, since the projection onto the brane of a null 
5D graviton momentum vector is in general a 4D timelike momentum vector. The massless mode arises when the 
projection is null, and corresponds to the 4D graviton, which is dominant at low energies. The massive Kaluza-Klein 
(KK) modes produce corrections to the gravitational potential in the weak-field static limit, which are 0(£^/r^) [5,6]. 

■ Table-top experiments currently impose an upper bound < 0.1 mm. 

■ On the brane, the negative As is offset by the positive brane tension A, and the effective cosmological constant on 
^) . the brane is 

A = i(A5 + «;2A), (1) 

U '. 

where = Sn/M^ and M4 ~ 10^^ GeV is the effective Planck scale on the brane. Fine-tuning can set A = 0. 
^ \ Because of the large extra dimension, the true fundamental gravity scale can be as low as ^ TeV in some brane-world 
scenarios, but in generalized RS2 models it is higher: 



5h ' /Vf2 
C3 ; Ml = -^, e< 0.1 mm ^ M5 > 10^ TeV , A > (100 GeV)^ . (2) 

At high energies {p 3> A) in the early universe, gravity becomes 5-dimensional and there are significant corrections to 
standard cosmological dynamics. 

The unperturbed cosmological brane-world is a Friedmann brane in a Schwarzschild-AdSs bulk [2] (where the ex- 
pansion of the universe can be interpreted as motion of the brane in the bulk) , satisfying the usual energy conservation 
equation but modified Friedmann equations. For a spatially flat universe without cosmological constant, these are 

n / p \ m 

A.4,a + „,(i + £)-2^i. (4) 

where w = p/ p and m is a constant, proportional to the mass of the bulk black hole. To avoid a naked singularity 
in the bulk, we take m > 0. The tidal. Coulomb effect of the 5-dimensional black hole on the brane is an effective 
radiative (cx a"'') term, the so-called "dark radiation". Nucleosynthesis is sensitive to the inclusion of additional 
relativistic energies not thermally coupled to the radiation plasma, and so is the cosmic microwave background. This 
places limits on the amount of dark radiation compatible with observations [7-9] 
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SO that p* -C Prad always, and the dark radiation is negligible by the matter era. The quadratic correction terms p^/X 
are dominant at high energies, p » A, and negligible at low energies. Since A » 1 MeV, the quadratic effects are 

negligible by nucleosynthesis. At late times, the standard 4-dimcnsional evolution is therefore recovered. 
In the absence of dark radiation, the exact solution (with ti; = = m) is [1] 

a = const [t{t + , tx = ^ < lO'^ sec . (6) 

VttX 

In the radiation era, w = ^, the exact solution with dark radiation {in ^ 0) is [7] 

a = const [t(t + tx)]'^^^ , tx = — , v^^4 

At late times, t'3>f\, wc recover the standard solution a oc ^^/^(w+i) ^ whereas at early times (high energies), t <^tx, 
the evolution is very different: a oc t^/^('^+^\ 

The high-energy brane-world correction provides increased Hubble damping in inflation; instead of iJ cx VV as in 
general relativity, the modified Friedmann equation (3) shows that H oc V when V ^ X. This means in particular 
that slow-roll inflation is possible even for potentials V{<f) that would be too steep in standard cosmology [4,10]. In 
addition, large-scale scalar perturbations generated by high-energy slow-roll inflation have an enhanced amplitude As 
compared with the standard general relativity amplitude (As)gr [4]: 

The background dynamics of brane-world cosmology are simple because they are effectively 4-dimensional. Cosmo- 
logical perturbations introduce truly 5-dimensional degrees of freedom, so that the 5-dimensional bulk perturbation 
equations must be solved in order to solve for perturbations on the brane. These 5-dimensional equations are partial 
differential equations for the 3-dimensional Fourier modes, subject to complicated initial and boundary conditions. 
However, on large scales, if one can neglect gradient terms, the density or curvature perturbations can be determined 
without knowing the 5-dimcnsional solutions [4,11,8]. The qualitative 5-dimcnsional effects on large-scale density 
perturbations during inflation and the transition to radiation-domination were investigated via a toy model in [12]. 
Here we generalise the results of [12] to include the dark radiation term m/a^ in the background. The presence of 
this term in the background can have a significant effect on the perturbations. 



II. FIELD EQUATIONS 

The effective field equations on the brane are derived from the 5D field equations in the bulk, ^^''Gab = — A5 '^^^gAB, 
by projecting the 5D curvature (using the Gauss-Codazzi equations), and then imposing the Darmois-Israel junction 
conditions at the brane (with Z2-symmetry) [3]: 

^2 

Gab = -^9ab + K^Tab + 6y5„6 - Sab ■ (9) 

Here Sab ~ (Tab)^ is the high-energy correction term, which is negligible for p <^ X, while Sab is the projection 

of the bulk Weyl tensor, encoding corrections from 5D graviton (KK) effects (and giving the dark radiation in the 
background). From the brane-observer viewpoint, the energy- momentum corrections in Sab are local, whereas the KK 
corrections in £ab are nonlocal [11,13,14], since they incorporate 5D gravity wave modes. These nonlocal corrections 
cannot be determined purely from data on the brane, and so the effective field equations are not a closed system. 
One needs to supplement them by 5D equations governing Sab, which are obtained from the 5D Einstein and Bianchi 
equations [3]. 

The trace free Sab contributes an effective KK energy density p*, pressure p*/3, momentum density q* and 
anisotropic stress tt*^ on the brane, 

\Sab = P* (uaUb + ^habj + Q^Ub + q^Ua + <(, , (10) 
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where u° is a physically determined 4- velocity on the brane and hab = ffafc + WaUb projects into the comoving rest space 
at each event. The brane "feels" the bulk gravitational field through these terms. In the background, g* = = tt*;,, 
since only the dark radiation term p* (oc ma^^) is compatible with Friedmann symmetry. The KK momentum density 
defines a velocity Va of the Weyl "fiuid" relative to by g* = p*Va- 

For a perfect fluid or scalar field, we choose w° as the frame in which there is no energy flux. Then the effective 
total energy density and pressure are 



eff . 



P 



P 

eflF 



P 1 + 
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Pi 
P 



- + (^ + 2-)|^ + f^. 



(11) 
(12) 



Energy-momentum conservation, V'Tob = 0, holds on the brane. Together with the 4D Bianchi identity, this implies 
that WSab = ViSoh /A, which shows qualitatively how 1+3 spacetime variations in the matter-radiation on the 
brane can source KK modes. These KK "conservation equations" in linearized form are 



P 



4 



{P + P) 
A 



(13) 
(14) 



where O is the volume expansion rate (= 3if in the background), Aa is the 4-acceleration, and is the covariant 

spatial derivative in the comoving rest space. Spatial inhomogencity (DaP 7^ 0) is a source for KK modes. Qualitatively 
and geometrically this can be understood as follows [11,13]; the non-uniform 5D gravitational field generated by 
inhomogeneous 4D matter-radiation contributes to the 5D Weyl tensor, which nonlocally "backreacts" on the brane 
via its projection £ab- Equation (14) shows that the source term is suppressed at low energies, and during quasi-de 
Sitter infiation on the brane. 

Equations (13) and (14) arc propagation equations for p* and g*. There is no propagation equation on the brane 
for TT*;,, so that one cannot determine the KK modes purely from data on the brane. In the 1+3-covariant description 
of braneworld perturbations [11,15-17], the KK anisotropic stress tt*;, is isolated as the term that must be determined 
from 5D equations. Once tt*^ is determined in this way, the 1+3 perturbation equations on the brane form a closed 
system. The KK terms act as source terms modifying the standard general relativity perturbation equations, together 
with the local high-energy corrections. 



III. DENSITY PERTURBATION EQUATIONS 

We define density and expansion (velocity) perturbation scalars, as in general relativity [18], 



(15) 



and then define dimensionless KK perturbation scalars [11], 



u = -dV* , Q = -D^g* , n = -B\* 



(16) 



where the scalar potentials q* and tt* are defined by g* = Dag*, tt*;, = (DaD(, — i/iabD^)7r*. The dark radiation 
fluctuation U is present even if there is no dark radiation in the background (p* = 0). It leads to an isocurvature 
(non-adiabatic) mode, even when the matter perturbations are assumed adiabatic [12]. We define the total effective 
dimensionless entropy S'^^ via the non-adiabatic part of the effective pressure as [18] 



where c^g = jf^ /jf^. Then 



[icl - 1 + (2 + 3w + 3c^) p/\] 

[(1 +w){l + p/X) + 4pV3p][3u; + 3(1 + 2u;)p/2A + p* / p] 



Ap* 
37 



A -(1 +«;)[/ 



(17) 



(18) 
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where = jyj p. If p* = in the background, then U is an isocurvature mode: S^'^ oc f/. If p* 7^ in the background, 
then the weighted difference between U and A determines the isocurvature mode: S'^'^ oc (4p*/3p)A — (1 + w)U . At 
very high energies, p 3> A, the entropy is suppressed by the factor A/p. At low energies in the radiation era (c^ = |) 
it is suppressed by the factor p/A. 

The covariant density perturbation equations on the brane [11] reduce to 



A = 3u;i7A - (1 



■ w)Z 

r.2 



1 + W 

U={3w- l)HU - 



D^A 



1 + 1 



K pU - -K p 



Q={3w- l)HQ 



3a 3 3a 



1 

4p* 
4c2 



l + w 
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Z - oD^Q , 
--3(l+«;)- 



l + w 



(19) 
(20) 

(21) 

(22) 



where / = u°'daf, and we have corrected the minor errors in the equations given in [12]. The KK anisotropic stress 
term E occurs only via its Laplacian, D^II. If we can neglect this term on large scales, then the system of density 
perturbation equations closes on super-Hubble scales [11]. An equivalent statement applies to the large-scale curvature 
perturbations [8]. We note that in the low-energy approximation scheme used by Koyama [19], 

HocD^C/, (23) 

so that in particular 11 ^ on large scales. 

On large scales in the braneworld, neglecting gradient terms, density perturbations are then described by the set 
{A, Z, U, Q}, where Q decouples from the other variables. KK effects introduce two new isocurvature modes on large 
scales, associated with U and Q [12,15]. 

We can find a first integral of the large-scale system by following the same approach used in general relativity [18]. 
A covariant local curvature perturbation is defined by the gradient of R^, the 3-Ricci scalar for u"-observers. This 
scalar is given by the Codazzi equation on the brane [14], 



R± = _^e^ + 2K^p-«. 

o 

Then the dimensionless curvature perturbation scalar is defined by C = a'^D^R^. This gives 

C = -Aa^BZ + iK^a^p (l + ^) A + In^a^pU . 

Using the large-scale equations (19)-(21), we find that C is locally conserved (i.e., constant along fundamental 
worldlines) on large scales. 



(24) 



(25) 



C = Co, Co = . 



(26) 



(The minor errors in [12] led to the erroneous conclusion that C was not in general locally conserved on the brane.) 

Thus C leads to a first integral for the large-scale system. We also replace A by the dimensionless perturbation 
variable 



$ = K^a^pA . 



(27) 



In general relativity, <I> is the covariant analogue of the Bardeen metric perturbation variable $h [18]. In the brane- 
world, high-energy and KK effects mean that the analogue of is a complicated generalisation of Eq. (27) (see [15]). 
The key problem in trying to use the generalised variable is that it contains 11, which remains undetermined. (This 

is the same reason that the large-scale Sachs- Wolfe effect on the brane cannot be determined from the density or 
curvature perturbations [8].) In any case, $ as defined above leads to useful simplification of the perturbation 
equations. 

The coupled system of Eqs. (19)-(21), which is closed on large scales (neglecting gradient terms), reduces via 
Eqs. (26) and (27) to the system 



$ = -if 



U = -H 



l + {l+w) 



l-3w 
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(l + w) 



3H^ 



U 



2p* 
Sa^Hp 



2 4 2 

a up 
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u + 





r P* 1 




Za^Hp 



Co ■ 



Co. 



(28) 
(29) 



(1 + w)k'^p 

This closed system for U generalises the equations given in [12] to the case p* ^ 0. Once $ and U are solved for 
Q is directly determined by Eq. (22) with D^II set to zero. 
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IV. SOLUTIONS 



The closed system for $ and U may be rewritten on using the Friedmann equation (3) and changing the time 
variable to the number of e-folds, t ^ N = ln(a/ao). This yields 



U' = - 



1 + 



3(1 + «;)(! + /)/ A) 



+ 



l-3«; 

P* 



2{l + p/2X + p*/p) 
2p*/p 



3(1 + 



2 2 

a K p 



il+p/2X + p*/p) 
1 1 



2{l + p/2X + p*/p) 

(1 + p/A) 



U- 



3(1 + u>) 



4(l + p/2A + p*/p) 



Co ) 



u- 









2a 



(1 + p/2A + p*/p) (1 + w) j a^K^p 



p {l + p/2\ + p*/p) a'^K^p 
where /' = df/dN and 



Co : 



a K p = P exp — / (1 + 3w)dN , (3 = a^n po . 



(30) 



(31) 



(32) 



Here /3 is dimensionless and can be given an arbitrary positive value by suitable choice of oq. 

For comparison, in the general relativity case (p/A = = p* = U) these two equations reduce to the single equation 



gr 



If w is constant, then the non-decaying attractor solution is 

_ 3(1 +w) 
^"^ 2(5 + 3w) 



Co. 



(33) 



(34) 



Deviations from general relativity arise from early-universe high-energy effects (when p ^ A) , and from the KK effects 
encoded in p* and U. At late times, the high-energy corrections become negligible. The background dark radiation 
decays like ordinary radiation {p* oc oc Prad), and is negligible at late times by Eq. (5). However, it can leave a 
signature that survives at late times, as we show below. 



A. No dcirk radiation in the background 



Equation (31) is much simpler in the case when there is no dark radiation in the background (p* = 0), i.e. when 
there is no bulk black hole. (Note also that the errors in [12] do not affect this case, so that the solutions given there 
are correct.) When p* = we get 

U = Uoexp- j{l-3w)dN . (35) 

Since p* = 0, it follows from Eq. (18) that U is an isocurvaturc mode. This isocurvaturc mode associated with dark 
radiation perturbations therefore decays for all equations of state less stiff than radiation, i.e. for w < ^. 
If w is constant, then Eq. (30) becomes 



(5-|-3«;) + (4-|-3«;)p/A 
2 + p/A 



3(l + «;) 



2(l + p/2A) 



-2N 



Uo 



3{l + w) 
4(l + p/2A) 



Co , 



(36) 



where Uq = (3Uq. The U contribution to $ is suppressed at low energies, even li w > ^ (when U grows). Thus the 
general relativity solution Eq. (34) for constant w is a low-energy attractor for p* = 0: 



low 



3(1 +w) 
2(5 + 3w) 



Co. 



(37) 



Thus there is no surviving imprint in $ of the dark radiation perturbations (the isocurvaturc mode) when p* = in 
the background. 

At high energies/ early times, Eq. (36) becomes 
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$'~ -(4 + 3w)$ 



-[3(1 + to)] e(i+3-)^t^o + - 

Po Po 



-(1 + w) 



(38) 



where A ^ po- Thus in the early universe, the contribution of the entropy mode U is sensitive to the equation of 
state. This contribution decays for to < — i, including the case of inflation. For w > Eq. (38) shows that the 
non-decaying solution is 



$high ~ — 

Po 



3{l + w) 



2(7 + 6m;) 
so that in the very high energy limit, 



Co 



^3{l+w)N _ ^ 

Po 



3/3(1 + w) 
(5 + 6m;) 



Uo 



3 A ,^ 
$high -> - — (1+m;) 
^ Po 



Co 



2C/n 



7 + 6m; 5 + Qw 



(39) 



(40) 



Hence, for constant w (> —5), high-energy effects initially suppress $ (by the factor X/po), before it starts to grow. 
This growth leads to the late-time general relativity value, Eq. (37). By contrast, in general relativity, $ remains 
constant when w is constant. 

The most relevant case is m; = |. If inflation takes place on the brane, reheating will initiate a radiation era. Provided 
the reheat temperature is high enough, the radiation era will start at high energies. Large-scale perturbations $ will 
be amplified according to Eqs. (40) and (37) with m; 
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N 

FIG. I. The evolution of $ for a mode that is well beyond the Hubble radius at iV = 0, about 50 e- folds before inflation 

ends, and remains supcr-Hubblc through the radiation era. A smooth transition from inflation to radiation is modeUcd by 
w = |[(| — e) tanh(A'^ — 50) — (f — e)], where e is a small positive parameter, chosen as e = 1/15 in the plot. (The parameter 
/? in Eq. (32) is set to I). Labels on the curves give the value of po/A; the general relativity solution is the dashed curve. 



For slow-roll inflation, w is not constant, but close to -I and slowly increasing, i.e., I + w ^ e, with < e <^ I and 
iJ-i|e| = |e'| < 1. Then Eq. (38) leads to 

^high-^eAcoe^^^, (41) 
^ Po 

so that $ has a slowly growing mode during high-energy slow-roll inflation. This is different from general relativity, 

where $ is constant during slow-roll inflation. Thus more amplification of $ can be achieved than in general relativity, 
consistent with Eq. (8). The evolution of $ is illustrated for a toy model of infiation-to-radiation in Fig. 1 (from [12]). 
The early (growing) and late time (constant) attractor solutions are seen explicitly in the plots. 



6 



B. Nonzero dark radiation in the background 



The presence of dark radiation in the background has interesting quaUtative impUcations, even though the quanti- 
tative constraints from observations on the amount of dark radiation are quite severe. The general relativity solution, 
Eq. (34) is no longer the late-time attractor when p* ^ 0. 

The dark radiation factor, for constant w, behaves as 



PL 
P 



ae 



{Zw-1)N 



a 



Po 



(42) 



where a is a dimensionless constant. For equations of state less stiff than radiation {w < i), the dark radiation factor 
1 

3' 



rcdshifts away. If w > i, then this factor grows with expansion. In the radiation era, w = |, the dark radiation factor 



is constant. 

If there was a period of inflation on the brane before the radiation era, then dark radiation would be redshifted to 

negligible levels. However, the reheating era that creates the radiation, could at the same time create dark radiation 
via 5D graviton emission in high-energy interactions [20]. The fraction of dark radiation is limited by observational 
constraints, Eq. (5), so that 



a < 0.03. 



(43) 



Wc can therefore perform a pcrturbative solution of Eqs. (30) and (31), up to first order in a. The zero-order solutions 
are given in the previous subsection. 

We find that the non-decaying solutions for radiation in the low-energy regime are 
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Uq + a— e 



2N 



(44) 

(45) 



It follows that the first-ordor correction to U (where U = U^^^ + U^^^) behaves like the large-scale density perturbations 
for radiation in general relativity: 



(46) 



In the high-energy regime, the non-decaying solutions are 
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C/high « C/q + a 
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8N 



(47) 

(48) 



In the very high energy limit. 



high 



high 



po 
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Co_ 
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539 



(49) 
(50) 



Comparing Eq. (49) with Eq. (40) shows that the high energy limit with dark radiation differs only negligibly from 
the no-dark-radiation limit. 
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FIG. 2. The evolution of $ on very large scales in the radiation era, with po/X — 10^ and dark radiation factor a — 0.03. 
Initial conditions are given by Co = 10~^ = Uo, and <1> is initially set to its high-energy attractor, Eq. (49). The solid curve is 
the numerical integration. The dashed curve is the high-energy approximation, Eq. (47). The dotted curve is the low-energy 
approximation, Eq. (44), and the dot-dashed curve is the low-energy approximation without dark radiation {a — 0). 

This analysis shows that $ is initially suppressed, then begins to grow, as in the no-dark-radiation case, eventually 
reaching an attractor which is less than the no-dark-radiation attractor. The role of dark radiation is to suppress, 
by the small fractional amount a, the final value of $, given in Eq. (44). The dark radiation perturbation U starts 
from a constant value at very high energies, and then grows at late times like the density perturbations in radiation. 
Figures 2 and 3 confirm the qualitative analysis. 
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FIG. 3. The evolution of U, with the same conditions as in Fig. 2. The dashed curve is the high-energy approximation, 
Eq. (48), and the dotted curve is the low-energy approximation, Eq. (45). 

The effective total entropy, Eq. (18), has the following behaviour, to 0{a), in the radiation era: 
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(1 - 2a)Uo . 



(51) 
(52) 
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At high energies \S'^^\ is suppressed by the factor X/po, and the a-correction is negUgible. At low energies, the large 
factor po/X is overwhelmed by the redshift factor e^*^. 

Without dark radiation in the background, U in the radiation era remains constant. When there is dark radiation 
in the background, U grows. The growth of U leads to a decrease in the final value of 

An extreme case of loss of power in $ occurs when the dark radiation rcdshifts more slowly than the matter, so that 
it eventually dominates over the matter. This happens when the matter equation of state is stiffer than radiation, 
w > For example, for stiff matter, w — 1, we have 

— cx , (53) 
P 

by Eq. (42). As a result, the decrease in <i> becomes so pronounced that $ is driven to at late times. By Eqs. (30) 
and (31), we find that 

[/lowcxe^^, $iow(xe-2^, (54) 
at late times. The suppression of $ is confirmed by numerical integration, as shown in Fig. 4. 




01 2345678 

N 

FIG. 4. The evolution of <1> on large scales in a stiff (kinetic) era (with parameters and initial conditions as in Figs. 2 and 3). 

This example could be relevant to the steep inflation scenario on the braneworld. High-energy effects mean that 
steep potentials can drive inflation on the brane when they would not give inflation in general relativity. Typically, 
this inflation ends in a kinetic regime with w ~ I, and with particle production being gravitational [10]. If dark 
radiation can be produced by high-energy interactions, then its presence serves to suppress $ on large-scales, as in 
Fig. 4. 

V. CONCLUSION 

We have investigated the effects of dark radiation in the background on large-scale density perturbations (neglecting 
gradient terms) in a Randall-Sundrum 2 type brane-world. Dark radiation significantly complicates the perturbation 
equations. Since it is limited to be a small fraction of the radiation energy density, its effects can be qualitatively 
analysed via a perturbative approach. This approach, confirmed by our numerical simulations, shows that in the 
radiation era, the large-scale density perturbations are suppressed at late times by a small amount. At the same time, 
the large-scale perturbations in the dark radiation itself grow at late times. By contrast these perturbations have 
constant amplitude in a universe with no dark radiation in the background. The suppression of density perturbations 
becomes strong when the matter has a stiff (kinetic-dominated) equation of state. 
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